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New class of symmetries for self-gravitating hydrodynamics 

equations 

Souichi MuRATA *•* and Kazuhiro Nozaki 
Department of Physics, Nagoya University, Chikusa-ku, Nagoya, 464-^602, Japan 



A method of calculating a new class of symmetries is presented for partial differential 
equations. The method give a new dynamical solution for an isothermal and cylindrically 
symmetric hydrodynamics equations under self-gravity. The solution describes a develop- 
^ |. ment of a gas cylinder from a motion-less state to states of expansion and oscillation. 



§1. Introduction 



Constructing exact solutions to a system of partial differential equations (PDEs) 
C/^ \ is an important problem in mathematical physics. One of the most powerful methods 

pH . to find particular solutions to a system of PDEs is based upon the analysis of its 

invariance with respect to the Lie group.^''^''^ If a system of PDEs is invariant 
under a Lie group of point transformation, we can construct special solutions, called 
similarity solutions, which are invariant under the group admitted by the system. 
This method provides a systematic and unified procedure in search of similarity 
solutions. 

fSJ ■ However, many systems of non- linear PDEs in mathematical physics admit only 

\ trivial symmetry groups, for example, translational symmetries with respect to time 

and space, rotational symmetries and scaling symmetries and so on. Invariant solu- 
Tj" ' tions under such trivial symmetry groups are also constructed by other naive methods 

\ based on physical considerations. If we find a method to derive a new non-trivial 

symmetry, a new type of self-similarity solutions may be constructed by means of 
the symmetry. 

In this paper, we present a method to obtain a new class of symmetry groups 
' admitted by a system of PDEs. A key of our treatment is to carry out symmetry 

^ ■ group analysis on particular solutions of physical interest. To illustrate our method, 

^ \ we apply the method to the equations of cylindrically symmetric isothermal hydrody- 

namics under self-gravitj^ . A Similarity solution for the self-gravitating isothermal 
gaseous cylinder without rotation has been discussed by Inutsuka and Miyama^' . In 
their study, they gave a self-similarity solution for a collapsing or expanding isother- 
mal filament whose density follows the Emden-type equation with an inertial term. 
Here, we study dynamics of the filament in the presence of a rotational flow around 
the axis of the cylinder by mean of our method of Lie group and find a new class of 
symmetries and a new quasi-self similarity solution with a rotational flow. 

The plan of the paper is as follow: in section 2, we explain our method to find 
a new class of symmetries to a system of PDEs. In section 3, we apply our method 
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to investigation of a quasi-similarity solution in the presence of a rotational flow. 

§2. Mathematical description of symmetry group analysis 

A system of n-th order PDEs in p independent and q dependent variables is 
given as a system of equations 

Z\^(x,u(")) = 0, u = l...l, (2-1) 

involving x = (x^, . . . , x^), u = {u^, . . . , u*?) and the derivatives of u with respect to 
X up to n, where u^") represents all the derivatives of u of all orders from to n. 

We consider a one-parameter Lie group of infinitesimal transformations acting 
on the independent and dependent variables of the system 

5^ = x* + er(x,u)+0(e2), i = l,...,p, (2-2) 
= n^' + e<JJ(x,u) +0(e2), j = l,...,q. (2-3) 

The infinitesimal generator V associated with the above group of transformations 
can be written as 

V = r (x, u)d,. + '^■'■(x, u)a„. . (2-4) 

i j 

The Invariance of the system (I2-1() under the infinitesimal transformations leads to 
the invariance condition 



Pr(")V[Z\^(x,u("))] = 0, iy = l,...l, whenever Z\/,(x, u^")) = 0, (2-5) 

where Pr^"^^ is the n-th order prolongation of the infinitesimal generator. 

Since the coefficients of the infinitesimal generator do not include derivatives of 
u, we can separate (|2-5() with respect to derivatives and solve the resulting overde- 
terminied system of linear homogeneous PDEs known as the determining equations. 
Our purpose for calculating symmetries is to obtain similarity solutions which are 
invariant under the symmetries. Applying the criterion of invariance, we get the 
invariant conditions called the functional self-similarity conditions 



u 



(2-6) 



where is the solution of the system. In order to solve the conditions H2-6p . we 
have to integrate the Lie equation corresponding to V 
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dx^ dx"^ dv} dvP' 

— • • • — — — ■■■ — ( • J 

Then, Eq H2-6() and solutions of 1)2 -7^ yield an invariant (self-similarity) solution 
under the symmetry. 

The crucial point is how to find a new non-trivial symmetry of the given equa- 
tions. Our approach is to restrict solutions on some special class, for example, in 
the following application we consider the Emden-type solution where the pressure is 
nearly balanced with the gravity. The restriction is given as an additional system of 
PDEs in general 



Z(x,u(")) = S(x,u), (2-8) 

where /^(x, u) is a function of x and u and will be given in calculating the determining 
equations. Then, the number of the determining equations is reduced owing to 
addition of the auxiliary equation H2-8p so that a new class of symmetries may emerge. 
To demonstrate how we can construct a non-trivial symmetry and a new solution, 
we apply the our method to the self-gravitating hydrodynamics equations. 

§3. Self-gravitating hydrodynamics for polytropic gas cylinder 

The self-gravitating hydrodynamics equations of a polytropic gas is often intro- 
duced in studies of the different stages of stellar evolutiorP . It often happens that 
the dynamics of a self-gravitating system exhibits self-similar behavior. A new simi- 
larity solution obtained either analytically or numerically, although constituting only 
a limited part of the general physical solution, may be very useful in understanding 
the basic dynamical behavior of the system. 

An isothermal polytropic gas cylinder is described by the following normalized 
fluid equation. 

0, (3-1) 
0, (3-2) 

0, (3-3) 

P, (3-4) 

where the pressure p , density p , radial velocity u , angular velocity v , gravitational 
potential a and radius r are normalized by po, po,UQ = vq = {"fPo/ PoV^'^ ,lPo/ Po^i^o, 
respectively; G is the gravitational constant and p = p"'. The acoustic time scale 
vq/uq, where tq = jpo/{4TrGpQ) , is set to equal to the time scale of free fall 
(AnGpo)-'/'. 
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Since we are interested in the Emden-type solution where the pressure is nearly 
balanced with self-gravity, we impose the following additional condition 



<7.r H Pr = -S{r,t), 

Pr 

or 

ut + uur = S{r, t), (3'5) 

r 

where S{r, t) is a function of r and t. 

Let us seek a Lie symmetry admitted by the system of Eas. ()c{-l|) . (|3-2|) . ()3-3|) . 
(|3-4|) and (|3-5j) . The function S{r,t) should be given by the determining equations 
for the symmetry. Applying Pr^^^V to Eqs. (12111), (E31), (EU and we 

obtain a dozen of determining equations for the symmetry. After some straightfoward 
calculations of the determining equations, we obtain the following symmetry. 



V = T{t)dt + -nrdr - Ttpdp + -^{-nu + Tttr)du - -nvd^ - -Tt(t)d^, (3-6) 

TutT - 2aTt = 0. (3-7) 
5(r,t) = ^r, (3-8) 

where (p = ar and a is an arbitrary constant. The invariance under the symmetry 
group (|3-6|) leads to the conditions of functional self-similarity on Eqs. (|3-1|) . (|3-2|) . 
(EH, (Eil 



-Tt/> - -^Ttrpr + t( + p(ur + - ) ) = 0, (3-9) 



1 1 / v"^ 1 

-(-TfU + Titr) - -nrUr +T{ur h -pr + ) = 0, (3-10) 

2 2 \ r p 

11 / 1 , , , 

--Ttf - -nrVr + TluVr + ) = 0, (3-11) 
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Let us choose the radial velocity as 



-lTt(l)-lTtr(p--(j)]-T(l)t = 0,. (3-12) 



u = ^r, (3-13) 

so that Eqs. ()3-9() and (|3-11|) are satisfied automatically. Then, Eqs ()3-5() . ()3-8|) . 
(nrm give 



2 4 6 8 10 12 

radius 

Fig. 1. density profile for a = 0, a = 0.05 and a = —0.05 with the boundary condition R{0) = 1, 

Ry{0)^0 



2t , [2t 2n 



(3-14) 



It is easy to see that Eq. (|3-7)) takes the following integrable form 



;;„, + 3A/,.-2»^.0. wtae /W^V7. (3-15) 
which can be integrated once and yields 

/« = -^, Fee A_alog(/), (S-IG) 

where b is an integral constant. 

Solving the Lie equation corresponding to the infinitesimal generator (|3-6|) . we 
get a new solution 



y = y' P=^r^ u = fty, V =0^2' '?^ = ~^' 

where the integral constant b should be positive. According to Eas. (|3-10|) and (|3-12j) . 
the function R{y) and ^{y) are expressed as the solution of the Emden-type equation 

54 + -3^ + exp(0) + 2a = O, cP = --^ - ay, e = log{R), (3-18) 
ay^ y ay ay 

where a term 2a comes from the inertia. 
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Fig. 2. Radial velocity for b = 0.1, b = 0.5 and b = 1.0 . 



For a = b = 0, the symmetry (|3-6|) becomes a naive Lie point symmetry admitted 
by Eqs. (|H-H) . (|:-{-2j) . (jSISI), l|»'{-4|) and the solution (|;-i-17|) coincides with a trivial self- 
similarity solution. For a ^ and 6 = 0, the solution ()3-17|) corresponds to the 
Inutsuka-Miyama solution. For other values of a and b ,the solutiion represents a 
new quasi-similarity solution. 

In Fig.l, typical solutions of the Emden-type equation H3-18() are plotted in the 
case a = 0, a = —0.05 and a = 0.05. A constant y , i.e., y = r/f{t) = y^ gives a 
radially moving frame, through which the gas dose not flow since u\y=y^ = dr / dt\y=y^. 
So, we can put the boundary of the density (the surface of the gas cylinder) on 
y = yo ■ The dynamics of the boundary depends on the evolution of f{t) only, which 
is determined by the potential V{f) given in Ea. (|3-16() 

Let us consider an initial value problem for the Eq. (|3-15|) such that f(t)\t=o = 1 
in order that y\t=o = r. The evolution of f{t) is classified into three cases. 

First of all, we consider the case a > and b > 0. In this case, we can take a 
motionless state initially so that ft\t=o = and ftt\t=o = a + b. Then, the surface is 
accelerated due to the positive inertia force. Radial and angular components of the 
velocity at the surface are depicted for a = 1 and b = 0.1,0.5, 1.0 in Fig.2 and 3. 

Next we consider a radially decelerating case without a rotational flow (i.e a < 
and 6 = 0), where we take ft\t=o = vo{> 0) as an initial condition so that the initial 
velocity of the boundary is positive. Then an initially expanding surface gradually 
decelerates due to the negative inertia term: Ea. H3-16|) can be integrated once and 
reads 




(3-19) 



A reflection point fr appears at 



Fig. 3. Angular velocity for b = 0.1, b = 0.5 and 6 = 1.0 . 
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Fig. 4. Radial velocity of the boundary for a = —0.5. 



/ = /. = exp^^ 

The surface stops at the reflection point and then the gas cylinder began to shrinks. 
Reflection time tr is given as 



i: 



V2(^o-Hlog(/)) 



In this reflection case, the radial velocity of boundary is depicted for a = —0.5 and 
i^o = 1 in Fig-4, where t,. = 4.06. 

Finally, we consider an oscillation case a < and 6 > 0. The potential V{f ) has 
a concave form and has an equilibrium point at / = /e = y6/|a| . Let us estimate a 
frequency around the equilibrium point, we expand / for small ^ 



Fig. 5. Radius / for a = 



-1.0, b = 2.0 and t;o = 




where ^ ^ 1 



(3-20) 



Substituting this expansion into Ea. H3-19() . we have 



2\dt 



) 



2 




eo-b 



(3.21) 



where the suffixation e of the each term denotes one at the equihbrium point and 
eo = Eo — V{y/b/\a\). Therefore, the displacement ^ is given as a harmonic oscillation 
^ = Asin^Lot) where 



In Fig. 5, 6 and 7, we show numerically calculated oscillating patterns /, the radial 
velocity and the angular velocity for a = —1.0, b = 2.0 and vq = 0.0. For this values 
of the parameters, the period is 7.23, which is larger than Itt/lo. 



We present a new method to the Lie symmetry analysis for non-linear PDEs. 
The key of our method is to calculate the Lie symmetry on special solutions. In 
mathematical physics, we often seek special solutions describing some physical phe- 
nomena. Characterizing a special class of solutions by an additional system, we find 
a broader class of symmetries. 




§4. summary 
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Fig. 6. Radial velocity for a = —1.0, 6 = 2.0 and vo = ■ 
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Fig. 7. Angular velocity for a = —1.0, 6 = 2.0 and vo = . 

To illustrate our method, we consider the Emden-type solution in self-gravitating 
hydrodynamics for an isothermal gas cylinder in the presence of rotating flow and 
find a new class quasi-similarity solution. 

The solution is classified by two parameters a and b, namely, the accelerating 
case (a > and 6 > 0), the decelerating case (a < and 6 = 0) and the oscillating 
case (a < and 6 > 0). The frequency near the equilibrium point in the third case 
is given by the values of parameters a and b. 
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